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Abstract 

We state and prove a generalized adiabatic theorem for Markov chains and provide 
examples and applications related to Glauber dynamics of Ising model over Z d /nZ d . 
The theorems derived in this paper describe a type of adiabatic dynamics for 
norm preserving, time inhomogeneous Markov transformations, while quantum adia- 
batic theorems deal with £ 2 (C n ) norm preserving ones, i.e. gradually changing unitary 
dynamics in C n . 
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Introduction 

The long-term stability of time inhomogeneous Markov processes is an active area of 
research in the field of stochastic processes and their applications. See [G] and references 
therein. Adiabatic times, as introduced in [2J, is a way to quantify the stability for a 
certain class of time inhomogeneous Markov processes. In order for us to introduce 
the reader to the type of adiabatic results that we will be working with in this paper, 
let us first mention earlier results that were published in [2], thus postponing a more 
elaborate discussion of the matter until section [TJ 

Mixing time quantifies the time it takes for a Markov chain to reach a state that is 
close enough to its stationary distribution. For the discrete-time finite state case we 
will look at the evolution of the Markov chain through its probability transition matrix. 
See [3] for a systematized account of mixing time theory and examples. Let || ■ \\tv 
denote the total variation distance. 

Definition 1. Suppose P is a discrete-time finite Markov chain with a unique station- 
ary distribution it, i.e. ttP = tt. Given an e > 0, the mixing time t m i x (e) is defined 
as 

tmix(t) = inf {t '■ H^-P* — ttHtv < e, for all probability distributions ^} . 
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To define adiabatic time in its first and simplest form (that we will expand and gen- 
eralize a few pages down) we have to consider a time inhomogeneous Markov chain 
whose probability transition matrix evolves linearly from an initial probability transi- 
tion matrix PinMai to a final probability transition matrix Pfi na i ■ Namely, we consider 
two transition probability operators, Pmitiai and Pfinai-, on a finite state space Q, and 
we suppose there is only a unique stationary distribution ttj of Pfi na i- We let 



Ps (1 SjPinitial & Pfinai 



(1) 



We use (pQ) to define a time inhomogeneous Markov chain P± over [0, T] time interval. 
The adiabatic time quantifies how gradual the transition from Pinitial to Pfinai should 
be so that at time T, the distribution is e close to the stationary distribution 7rj of 

Pfinai ■ 

Definition 2. Given e > 0, a time T e is called the adiabatic time if it is the least T 
such that 

max u \\vPi_P2_ ■ ■ ■ Pt-i P] — \\tv < e 

rp rj-i rj-i J ' ' 

where the maximum is taken over all probability distributions v over £1. 

With these definitions one would naturally ask how adiabatic and mixing times 
compare. This will be especially relevant given the emergence of quantum adiabatic 
computation and some instances of using adiabatic algorithms to solve certain classical 
computation problems. See pQ and [5]. It can be speculated that there may be scenarios 
in which the adiabatic time is more convenient to compute than mixing times. If we 
find the relationship between the two, it will give us an understanding of the adiabatic 
transition (which is more prevalent in a context of physics) in terms of mixing times 
and vice versa. The following adiabatic theorem was proved in [2]. 



Theorem (Kovchegov 2009). Let t r 

adiabatic time 



O 



denote the mixing time for P final- Then the 



tjnix 



In subsection 11.21 we will give an example to show that t^ mix is the best bound for the 
adiabatic time in this setting. There f2 = {0, 1,2, ... ,n} and 
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Similar adiabatic results hold in the case of continuous-time Markov chains. There, 
the concept of an adiabatic time is defined within the same setting and a relationship 
with mixing time is shown. Let us state a continuous adiabatic result from [2], and 
then prove a more general statement of the theorem in the next section. 
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Once again we define mixing time as a measurement of the time it takes for a Markov 
chain to reach a state that is close enough to its stationary distribution. For the 
continuous-time, finite-state case we look at the evolution of the Markov chain through 
its probability transition matrix as a function over time. 

Definition 3. Suppose P(t) is a finite continuous-time Markov chain with a unique 
stationary distribution ir. Given an e > 0, the mixing time t m i x (e) is defined as 

t m i x {e) = inf {t : \\vP(t) — tt\\tv < e > f or a ^ probability distributions v} . 

To define an adiabatic time we have to look at the linear evolution of a generator for the 
initial probability transition matrix to a generator for the final probability transition 
matrix. Suppose Qinitiai and Q final are two bounded generators for continuous-time 
Markov processes on a finite state space f2, and irf is the unique stationary distribution 
for Q fi na i . Let us define a time inhomogeneous generator 

Q[s] = (1 - s) Qinitiat + sQ final (2) 

Given T > and < t\ < t% < T, let Pr(ti,t2) denote a matrix of transition 
probabilities of a Markov process generated by Q[-f] over the time interval [ti , ^2] • 
With this new generator we define the adiabatic time to be the smallest transition 
time T such that regardless of our starting distribution, the continuous-time Markov 
chain generated by Q[j^} arrives at a state close enough to our stationary distribution 

*f- 

Definition 4. Given e > 0, a time T t is called the adiabatic time if it is the least T 
such that 

max IIzaFWO, T) — 717 \\tv < e 

where the maximum is taken over all probability distributions v over £1. 

The above definition for continuous-time Markov chains is similar to the one in the 
discrete time setting. The corresponding adiabatic theorem for the continuous times 
case was proved in [2]. 

Theorem (Kovchegov 2009). Let t m % x denote the mixing time for Q final- Take A such 
that A > maxjgsi J2j : j& Qi™ tml and ^ > max i6 n Ylj-.j^i %™ 1 > where q fJ tml and q{j al 
are the rates in Qinitiai & n d Q final respectively. Then the adiabatic time 

j, < At m i x (e/2) 2 ^ 
e 

where 6 = t mix {e/2) + e/(4A). 

This is once again the best bound as can be shown through the corresponding example. 

In the next section we will state the adiabatic results for Markov chains that generalize 
the above mentioned theorems in [2] and provide examples of applications in statistical 
mechanics. Section [2] is dedicated to proofs. 
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1 Results and applications 

Here we extend the results from [2], and thus expand the range of problems that can 
be analyzed with these types of adiabatic theorems. One such problem that we will 
discuss in subsection 11.11 deals with adiabatic Glauber dynamics for Ising model. Now, 
in order to solve a larger class of problems, we redefine the adiabatic transition for 
both the discrete and continuous cases. 



We consider an adiabatic dynamics where transition probabilities change gradually 

from Pinitial = to P final = {W^} so that, for each pair of states i and j, 

the corresponding mutation of pij from p l P-} tial to p{ l ™ al is implemented differently and 
not always linearly. In the case of discrete time steps, this means defining 

„y W = (i - <M*)Kf a/ + hi&piT 1 > ( 3 ) 

where cpij : [0, 1] — >• [0, 1] are continuous functions such that </>ij(0) = and 4>ij(l) = 1 
for all locations (i,j). 

The above definition generalizes ([I]). If we suppose there is a unique stationary distri- 
bution TTf for Pfinai, then the Definition [2] of adiabatic time T e given in the previous 
section will hold for the adiabatic dynamics defined in ([3]). The new T e is related to 
mixing time via the following adiabatic theorem, that we will prove in section [2j 

Theorem 1 (Discrete Adiabatic Theorem). Let P± = {pij [y] } be an inhomogeneous 
discrete-time Markov chain over [0, T]. Let 4>(s) = minjj (pij(s) be the pointwise min- 
imum function of all of the 4>ij functions. If m > 1 is an integer such that (j) is m + 1 
times continuously differentiable in a neighborhood of 1, 

<^>( fe )(l) = for all integers k such that 1 < k < m 

and (p( m \l) / 0, then 



O 



m + l 

fc/2) 



6 m 



The above is, in fact, the best bound in the new setting as shown through the example 
given later. See subsection 11.21 

Now we extend the notion of adiabatic dynamics for the continuous-time Markov gen- 
erators as follows. We let 

qi>j [s] = (1 - Msfiqtf 01 + Ms)q{p l for all pairs i + j, (4) 

where once again faj : [0, 1] — > [0, 1] are continuous functions such that (f>i t j(0) = and 
(j>i } j(l) = 1 for all locations (i,j). Also, we let Q[s] denote the corresponding Markov 
operator. 

If there is a unique stationary distribution TTf for Qfi na i, then the Definition [4] of 
adiabatic time will apply for the extended adiabatic dynamics in ([4]), and the new 
T e can be again related to mixing time. 
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Theorem 2 (Continuous Adiabatic Theorem). Let Q [A] (t £ [0,T]) generate the 
inhomogeneous discrete-time Markov chain. Let <j)(s) = mini <j>i,j(s) be the pointwise 
minimum function of all of the 4>ij functions. Suppose m > 1 is an integer such that 
(j) is m + 1 times continuously differentiable in a neighborhood of 1, 

(f)( k \l) = for all integers k such that 1 < k < m 

and <j)( m \l) / 0. If we take A such that 

A > max V gff^ and A > max V g* na, J 

where q"V Ual and q{™ al are the rates in Qinitial and Q final respectively. Then 



O 




The reader can reference the proof of this theorem in section [2j Again this is the best 
bound in the new setting as can be shown through the same example. See subsection 

ma 



Now we check that the above continuous adiabatic theorem is scale invariant. For 
a positive M, we scale the initial and final generators to be jjQinitial an d jjQ final 
respectively. Then the adiabatic evolution is slowed down M times, and the new 

_L_ m + 1 

adiabatic time should be of order M [^1 m i m ™, (e/2) with the old t m i x and A taken 
before scaling. On the other hand the new mixing time will be Mt m i x , and the new A 
is as the rates are M times lower. Plugging the new parameters into the expression 
in the theorem, we obtain 



Me 



(Mt mix ) 



m + 1 



M 



m + 1 



confirming the theorem is invariant under time scaling. 



Let us revisit adiabatic theorems in physics and quantum mechanics. The reader can 
find a version of quantum adiabatic theorem in [3] and multiple other sources. 

The adiabatic results in physics consider a system that transitions from one state to 
another, while the energy function changes from an initial Hinitiai to H final- If the 
change in the energy function happens slowly enough, for the system that is initially 
at one of the equilibrium states (i.e. at an eigenvector of the initial energy function 
~Hi n itial)i the resulting state of the transition of the system will end up at a state that 
is e close to the corresponding eigenvector of the final energy function H final- That is, 
provided the change in the external conditions is gradual enough, the jth eigenstate of 
H-initiai is carried to an e proximity of the jth eigenstate of % final- 

Often the adiabatic results concern with one eigenstate, the ground state. Thinking 
of Schrodinger equation as an £ 2 (C n ) norm preserving linear dynamics, and a finite 
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Markov process as a natural description of an £ l (W+) norm preserving linear dynamics, 
the ground state of one would correspond to the stationary state of the other. It is 
important to mention that in addition to all above properties, the quantum adiabatic 
theorems often require the transition to be gradual enough for the state to be within 
an e proximity of the corresponding ground state at each time during the transition. 
Taking this into account, the complete analogue of quantum adiabatic theorem for 
£ 1 (R") would be the one in which the initial distribution is jiQ = n 'initial and 

||/xt-7r t || < e Vt€[0,T], 

where [it = Mo -Pi -Pi ■ ■ ■ Pt-i Pi is the distribution of the inhomogeneous Markov chain 

T T T 

at time t G [0, T], ^initial is the stationary distribution of Pinitial, and ir t is the stationary 
distribution P±. While we are currently working on proving the above mentioned 

T 

complete analogue in both discrete and continuous cases, the adiabatic results of this 
section are sufficiently strong for answering our questions concerning adiabatic Glauber 
dynamics as stated in the following subsection. 

1.1 Applications to Ising models with adiabatic Glauber 
dynamics 

Let us first state a version of the quantum adiabatic theorem. Given two Hamiltonians, 
^initial and % final, acting on a quantum system. Let 

%(s) = (1 - s)n in uial + sUfinal (5) 

Suppose the system evolves according to H(t/T) from time t = to time T . Then if T 
is large enough, the final state of the system will be close to the ground state of H final- 
They are e close in the £ 2 norm whenever T > where f3 is the least spectral gap 
of H(s) over all s 6 [0, 1], and C depends linearly on a square of the distance between 

T~Linitial and H final- 

Now, switching to canonical ensembles of statistical mechanics will land us in a Gibbs 
measure space with familiar probabilistic properties, i.e. the Markov property of sta- 
tistical independence. We consider a nearest-neighbor Ising model. There the spins 
can be of two types, -1 and +1. The spins interact only with nearest neighbors. A 
Hamiltonian determines the energy-value of the interactions of the configuration of 
spins. 

Here, for a microstate, we multiply its energy by the thermodynamic beta and call 
it the Hamiltonian of the microstate. In other words, letting x be a configuration of 
spins, the Hamiltonian we use in this paper will be defined as 

%{x) = -^Y i M i jx{i)x{j) 

where ft is the thermodynamic beta, i.e. its inverse is the temperature times Boltz- 
mann's constant, M = {Mjj} is a symmetric matrix and for locations i and j, Mij = 
if i is not a nearest neighbor to j and My = 1 if i is a nearest neighbor to j. 
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The Markov property of statistical independence is reflected through the local Hamil- 
tonian defined at every location j as follows 

n loc (x(j)) = -p *(0*0'), 

where i ~ j means i and j are nearest neighbors on the graph. 

In the original, non-adiabatic case, the Glauber dynamics is used to generate the 
following Gibbs distribution 

ir( x ) = 1 e -U(x) 

over all spin configurations x £ {— where S denotes all the sites of a graph, 
and Z(f3) is the normalization constant. Let us describe how the Glauber dynamics 
works in the case when each vertex of the connected graph is of the same degree. 
There for each location j we have an independent exponential clock with parameter 
one associated with it. When the clock rings, the spin x(j) of configuration x at the 
site j on the graph is reselected using the following probability 

P(x(j) = +1) = „„ ; ; rr -, , , rr = 2 - 2 tanh \n loc (x + (j)) } 

y KJJ ' e ~U loc {x-(j)) + e -H loc {x+(j)) \ y +\J))j 

where x + {i) = X-(i) = x(i) for i / j, x + (j) = +1 and = —1. Here 

P(x{j) = -1) = 1 - P{x{j) = +1) 

Also H loc '{x -(j)) = —l-L loc (x+(j)). Now we have a continuous-time Markov process, in 
which our states are the collection of the configurations of spins. 

Now, consider an adiabatic evolution of Hamiltonians as in ([5]). There at each time t, 

%{s) = (1 — sfH-initial + s1-Lfi na i, 

where s = A. The local Hamiltonians must therefore evolve accordingly, 

n s — \ L — s ) n initial ' sn final 

and the adiabatic Glauber dynamics is the one where when the clock rings, the spin 
x(j) is reselected with probabilities 

e -u i r{x + {j)) 

P s (x{j) = +1) = e _ H loc (x _ U)) + e _ H loc (x+U)) 

Here too, ^ oc (x_(j)) = -n{ oc (x + (j)). 

The stationary distribution of the Qfi na i -generated Markov process, i.e. Glauber dy- 
namics with Hinffial energy, is, for a configuration x, 

e —Wfina,l( x ) 

tt(x) 



Sail config. x' 
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1.1.1 Adiabatic Glauber dynamics on Z 2 jnl? 

Consider nonlinear adiabatic Glauber dynamics of an Ising model on a two-dimsnsional 
torus Z 2 /nZ 2 . There any two neighboring spin configurations x and y in {—1, +l} n dif- 

I x(vl\ if u v 

fer at only one site on the graph, say v € 1? jnl? . That is y(u) = < 

I 1 — x(u) if u = v 

The transition rates evolve according to the adiabatic Glauber dynamics rules, and the 
transition rates can be represented as 

q x ,y[s) = (1 - + ^MllT 1 

as in Here the functions 4> x ,y(s) for two neighbors x and y depend entirely on the 
spins around the discrepancy site v. Namely if all four neighbors of v are of the same 
spin (+1 or —1), then 

, , v cos/i(-4/3 2 ) • sinh(s(4,/3t - 4/3 2 )) 



' sinh(4Pi - 4/3 2 ) • cosh(-4fix + s(40i - 4/3 2 )) 

If it is three of one kind, and one of the other (i.e three +1 and one —1, or three —1 
and one +1) as illustrated below 

-1 

+1 - v - +1 
+1 

then 

cosh(-2p 2 ) ■ sinh(s(2/3i - 2/3 2 )) 



sinh{2(5i - 2fa) ■ cosh(-2^ + s(2/3 1 - 2/3 2 )) 



If there are two of each kind, any function works, as both, the initial and the final, 

init; 
ix,y 



local Hamiltonians produce the same transition rates q % ™y tal = ql]y al = 1/2. 



Now, since 



cosh(-A/3 2 ) ■ sinh(s(4/3i - 4/3 2 )) cosh(-2/3 2 ) ■ sinh{s(2p x - 2/3 2 )) 



sinh{A(3i - 4/3 2 ) • cos/t(-4/3i + s(4/3i - 4/3 2 )) ~ sinh{2f3i - 2f3 2 ) ■ cosh(-2fi x + s{2(3i - 2(3 2 )) 
for s € [0, 1], Theorem [2] implies the adiabatic time 

/ 2 \-i2 N 
log(n) + log ( - 

where C = ( 2 / 3l ~ 2 / 32 )[ c °*M 2 ^i^^)-fanfe(-2fe)] ^ jj er6) a t every vertex on the torus we 

attached a Poisson clock with rate one, and therefore we can take A = n 2 . Also m = 1 
in the theorem, and one can find the expression for t m i x in [3]. 
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1.1.2 Adiabatic Glauber dynamics on Z d /nZ d 

The adiabatic Glauber dynamics of an Ising model on a d-dimsnsional torus *L d jnL d 
solves similarly. There the minimum function of the Theorem [2] is same as in the 
case of d = 2 



(s) 



and the adiabatic time 



cosh(-2(3 2 ) ■ sinhjs^i - 2(3 2 )) 
sinh{2l3 x - 2(3 2 ) ■ cos/i(-2/3i + s(2f3 1 - 2(3 2 )) 



O C 



n 



log(n) + log 



where again C 



(2ft-2fe)[coffe(2/3i-2/3 2 )-tanfc(-2/3 2 )] 
{l-tanh(2(3 2 )] 2 ' 



Notice that the time scaling argument that followed the statement of Theorem [2] works 
here as well. That is, if we use one Poisson clock of rate one for all vertices, or equiv- 
alently place Poisson clocks of rate n~ d at every individual vertex, the new adiabatic 
time will be 



T' = 0\C 



n 



2,1 



log(ra) + log 



as A' = 1 here. 



1.2 The bound is optimal 

In this subsection we give examples that show that the ij^. order of adiabatic time 

m+l 

given in Kovchegov [2], and t ^ order for more general settings of this current paper 
are in fact optimal. We consider discrete parobability transition matrices 

••• \ 

••• 

1 '•• : 

■• '■• 
• • 1 
•• l) 

over n + 1 states, {0, 1,2,..., n}, and let the discrete-time adiabatic probability tran- 
sition matrix to be 

Ps — (1 S)Pinitial $P final 

as in [2j. Let Hf again denote the stationary distribution for Pfi na i. Here Hf = 
(0, . . . , 0, 1) and the mixing time t m i x (e) = n for any e £ (0, 1). 

Now, since iiPmiUai = P f° r an y probability distribution fi, we have the following 
inequality 



initial 











(° 


1 





( 1 


•• 


• °\ 
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and Pfi na i — 
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\uPt ■ ■ ■ Pt — Try 1 1 

rr, m " 



> 



3=0 



T\ 



T ) j! ■ T T ~i 



TV 



Observe that pPj in ? = Try for any < j < T — n. Therefore 
\\vPi ■ ■ ■ Pt ~ lT f \\TV > Yl i 1 '^ 



j=T-n+l 
T 

£ 

j=T-n 



T T -i (j - 1)! • T^-O'- 1 ) 

T! 



(T-n)\-T n 
T-n+1 T- 1 



Now, because T ~y +1 • • • < ( T 2 J for n > 2, we see that 



\vPi_ • • • Pt — 7Ty||ry > 1 



Thus e > jl^-Pj. • • • Pt — 7Tf\\TV > 1 ~~ e 



> 1 -e 



implies T > 



n? ti 



41og(l-e) 4e 



in / .r 

4e ' 



proving that the order of adiabatic time T e = O 



in [2] is optimal. 



1.2.1 Optimal bound for general settings 



The same example works in the more general setting considered in this paper. For the 
same P ini tiai and Pf ina i, let 



p i j[s] = (l- ( f> i j(s))ptj tial + (t> ijj (s)p. 



final 



as in (j3J). Suppose <t>i,j(s) = 4>(s) for all pairs of states (i, j), and suppose m > 1 is an 
integer such that eft is m + 1 times continuously differentiable in a neighborhood of 1, 

<^( fe )(l) = for all integers k such that 1 < k < m 

and ^ m \l) / 0. Then 



\vP± ■ ■ ■ Pt — TTf\\TV > 



Z=0 3=1+1 



TV 
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and therefore 

T-l T 

\\uP,...PT-n f \\ T v> £ Q.-W/T)) J! <Ki/T) 

i=T-n+l i=i+l 
T-l / T T 

= e n *oyT) - n ^^'/ t ) 

/=T-n+l \j=/+l j=l 
T 

=1- n 4>(j/T). 

j=T-n+l 



The minimum function <f>(x) = 1 + 

T 



0( m )(l)(x-l) r ' 



+ (Is - ll m+1 ) and 



vPy .. P ^, bv ^. n L '-'^'tf-" +0 ((i-i/r) 



m+l^ 



j=T-n+l 



> 



T m ■ ml 

1 _ E^-„ +1 lo g (l + t -" m ^^' T -^ + 0((l- J /T)^ 1 )) 

1 _ e (-l)y(l) . E n-l^ +0((n/r)m+ l ) 



as log(l + x) < x. 



It is a well know fact that 

n-1 



j=l j=0 



B, 



(k + 1) - j Vi 



O+i) -j 



where l?j is the jth Bernoulli number. 
Suppose e > \\vP±. • • • Pt — iTf\\Tv, then 



e « - log(l - e) > 



E 

i=o 



5, 



(m + 1) — j\j 



in 



n 



(m+l)-j 



(6) 



+ 0((n/T) 



Thus confirming that the order of adiabatic time T e = O 
optimal. 

Naturally, there is a similar example in the continuous case. There 



in Theorem [T] is 



Q 



initial 



( 
1 -1 

1 -1 



\l 



o \ 




-lj 



and Q final 
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0-110 

0-11 
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2 Proofs 

In this section we give formal proofs to both adiabatic theorems, Theorem [1] and 
Theorem [2j 



2.1 Proof of Theorem [T] 

Proof. We write 

Pitj [s] = (i - M8))p%™»> + {Ms) - <Ks))pVT' + Mp^ 

and define transition probability matrix P to be such that 

(i - <K*))p y = (i - MsVpff*** + (Ms) - Hs)) P Z inal) 



We will thus have that 



Observe that 



P s = (1 - <Ka))P + <f>(a)P final . 



VP1P2 •••Pt-iPi 



n wm 

j=N+l 



»NPjr» + £ 



where vn = vP^P^. ■ ■ ■ Pit, and £ is the rest of the terms, and both T and TV" are 

T T T 

natural numbers with N < T. 

By the triangle inequality, we have 



rp at 

max H^PiPi • • • Pr-i Pi — 7T/||tV < max \\ V P fined ~ ^fnTV • 



where < Sn < 1 



n 

j=N+l 



+ S N 



X\j=N+l Y\TJ 

Let set T — N = i m j x .(e/2), where e > is small. Then we have that 

T 



m ^\W P final -W/IItV • 



< e/2 



Setting 1 — nJ=Ar+i ^(40 — e /2 we obtain 



log (1- e/2) < l °S0U/T) 

j=N+l 

We plug in the approximation of the minimum function cj) around x = 1 

0(»»)(i)(x - ly 



(X =1 + 



m: 



+o(|x-ir +i ) 
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obtaining 



* <!-</*>>- E Ji t ( -"t l !"' r ^i'Hr')' 



j=7V+l \ 
Therefore 

_ log (1 _ tm > Lgg?> T f + o (y- N r 2 ) 

Observe that (-l) m+1 c/> {m) ( 1 ) > as : [0, 1] -> [0, 1] and <£(1) = 1. 

By©, E;=T (e/2) " 1 i m = Er=0(^^(T)^( e / 2 ) (m+1)_fc ' where 5 fc is the fcth 
Bernoulli number, and therefore 



e> log(1 c/2 w (-ir + v m) (D f M f . {c;2) ^-^ c ( {T - N)m+2 \ 



In order for the right hand side of the above equation to be — log(l — e/2) close to 

/ m±l \ 

I t m (e/2) \ 

zero, it is sufficient for T to be of order of O mix t . □ 



2.2 Proof of Theorem d 

Proof. Define Q to be a Markov generate with off-diagonal entries 

1 - <j>i,j(s) (initial) <t>i,j{s) - 4>(s) (final) 
qi >3 ~ 1 _ 0( s ) %3 + 1 _ <?i,j 

Then writing 

%j [s] = (1 - M (,))gJ7 tM ° + (M*) " ^))flg iBa0 + <P(s)Q% nal) 
would imply 



Observe that 



as 



Q[s) = {l-4>{s))Q + <j){s)Q final 



A > max > di ,• and A > max > % ,• 



s ^ (initial) , , ^_ (final) 

A > max > o> • and A > max > q- ■ 

~ ien ^ M ~ ien ^ hl 

Let Pfinaiit) = e t( ^f inal denote the transition probability matrix associated with the 
generator Qf ina u and let P Q = I + \Q and P X = I+ jQ fina i. 

The Po and Pi are discrete Markov chains. Conditioning on the number of arrivals 
within the [N, T] time interval 



\n=0 



uP T (0,T) = v N P T (N,T) =»n(J2 ^ T - N ))\ -KT-N) In 



Adiabatic times and applications 

where un = i/Pt(0,N) and 
n! 



(T - N) n 



N<s 1 <-<s n <T 



dsi • • • ds n 



i.e. the order statistics of the n. 

Therefore, combining the terms with P final-, we obtain 



vPt^T) = u n [J2 



final —X(T—N) 1,1 - / £l 

n! .L .L 



\n=0 



\n=0 



N JN 
1 



(?) 



c/si • • • ds n + £ 



n! final \ in 



4>{x)dx 



where £ denotes the rest of the terms. 



Take K > and define 



T= lj Ki <P(x)dx) t mix (e/2) 



and 



-l 



V =¥—^\ C o(x)dx I /„,,,(r/2) 



K \ Ik-l 

K 



Recall the approximation of the minimum function <fi around x = 1 

0M(l)(x - l) r 



<^(x) = 1 + 



m! 



+o(|x-ir +i ) 



and therefore 



where = (-l) m ^iy + 0(if _1 ). Thus we can write 



Vn=0 



n! 



final 



1 + 7(iT) 



T — N 



+ £ 



We see, using standard uniformization argument, that 



zaP t (0,T) = e~ 



A"t m ia;(e/2)"- \ 

; -^/w + £ 



vn=0 



n! 



3 ^ ( -ff "»+7(K) 



i/jvexp{Q/ irio j • i mix (e/2)} + £ 



Adiabatic times and applications 



Now, since (_l)m^M(i) < o, we have that, for any probability distribution u, 

\\uPr(Q,T) - 7T f \\ TV = \\vexp {Q final ■ t mix (e/2)} - ir f \\ TV • e A (^™)'— (e/2) + S N , 
where, by triangle inequality, 

' A"(t mix (e/2)) n 



< s N < i - e A ( l+ ^) 1 -*™«( e / 2 ) f jr ■ 

\n=0 



and, by definition of i m ix> 

||^exp{Q /w -^(e/2)} -tt/Htv • e A (^™K" {e/2) < e/2 
Taking K = c(A/e)™ t mix (e/2)™ with constant c >> (-l) m+1 f^+iy. > we obtain 

e > - log(l - e/2) > A f ^ ( 5U l 

and therefore 
Thus 



< S N < 1 - e K )* mia:(e/2) < e /2 



^ _ Kt m j x (e/2) _ 

1 | 7W " U 



( 


"A" 




e 



1 

m+1 



t m 7 x (f/2) 



□ 
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